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We report measurements of the n~p — » n°n differential 
cross sections at six momenta (104-143 MeV/c) and four an- 
gles (0-40 deg) by detection of 7-ray pairs from n° — > 77 
decays using the TRIUMF RMC spectrometer. This region 
exhibits a vanishing zero-degree cross section from destructive 
interference between s- and p-waves, thus yielding special 
sensitivity to pion-nucleon dynamics and isospin symmetry 
breaking. Our data and previous data do not agree, with 
important implications for earlier claims of large isospin vio- 
lating effects in low energy pion-nucleon interactions. 
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Studies of pion-nucleon dynamics have long played a 
central role in nuclear physics. For example, the 7rN sys- 
tem enables important tests of basic symmetries - no- 
tably isospin symmetry and chiral symmetry - that are 
fundamental to the low energy realization of the strong 
force. Interestingly, the possibility of gross violations of 
isospin symmetry in 7rN processes involving neutral pions 
was first suggested by Weinberg [1]. 

Of the eight possible 7rN-^7rN channels there are three 
experimentally accessible channels: 7r + p elastic scatter- 
ing, 7r~p elastic scattering and n~p — > ir°n charge ex- 
change (CEX). Assuming isospin symmetry, the three 
hadronic amplitudes for the three channels are deter- 
mined by two isospin amplitudes, thus implying a tri- 
angle relationship among the amplitudes fcEX = (/- — 
/+)/v2. Consequently, an important test of isospin sym- 
metry is possible by comparing the extracted amplitude, 
fcEXi from CEX data, with the predicted amplitude, 
f' CEX = (/_ — f + )/y/2, from elastic data. 

Special sensitivity to isospin violation in 7rN interac- 
tions is afforded by pion CEX at forward angles and 100- 
150 McV/c. At these kinematics the dominating s- and 
p-waves interfere destructively and yield a striking dip in 
the cross section [2] . This cancellation of the isospin con- 
serving amplitude thus amplifies the sensitivity to possi- 
ble isospin non-conserving amplitudes, thereby allowing 
a delicate test of symmetry breaking. 

A comparison between the scattering amplitudes pre- 
dicted from elastic data {fc EX ) anc ^ those extracted from 
CEX data (fcEx) was made by Gibbs et al. [4,5] and 
Matsinos et al. [6,7]. While different in details, both 
studies found evidence of an unexpectedly large isospin 
violation after accounting for the effects of 7r ± p Coulomb 



interactions and hadronic mass differences (we note this 
conclusion was questioned in Refs. [8,9]). In particu- 
lar, both Gibbs et al. and Matsinos et al. reported a 
—0.012 ±0.003 fm discrepancy between the forward scat- 
tering amplitudes that are extracted from the CEX data 
and those predicted from the elastic data. 

Subsequently, Piekarewicz [10] claimed that a symme- 
try breaking originating in the u,d-quark mass difference 
may explain these results. In Piekarewicz's calculation 
the u,d-quark mass difference induces inequalities among 
the 7r°p, 7r°n and 7r ± N couplings that are realized as a - 
0.015 fm correction to the triangle relationship. However 
Fettes and Meissner [11], using chiral perturbation the- 
ory, have concluded quark mass differences yield much 
smaller isospin breaking effects in the 7rN interaction. 

While the elastic dataset is extensive and precise, the 
CEX dataset is more limited and less accurate [12]. 
Most notably, the region of momenta 100-150 MeV/c 
and forward angles is dominated by a single experiment 
due to Fitzgerald et al. [2] (an additional experiment 
by Isenhower et al. [3] has reported preliminary data 
at 112 MeV/c). The Fitzgerald experiment used the 
LAMPF ir° spectrometer and measured cross sections at 
momenta 100.6-147.1 MeV/c and angles #<25 degrees. 
Unfortunately, both the absolute normalization, obtained 
from 12 C(7r,7rN) 11 C activation measurements, and the 
background subtraction, determined from 12 C(7r~,7r°) 
background measurements, have been questioned [3,13]. 

Given the importance of the cross section minimum as 
an isospin symmetry test we undertook a new measure- 
ment of CEX differential cross sections at momenta 104- 
143 MeV/c and angles 0-40 degrees. Our experiment was 
conducted using the RMC detector on the M9A beam- 
line at the TRIUMF cyclotron. The RMC detector is a 
large acceptance pair spectrometer that was developed 
for radiative muon capture studies [14]. 

The M9A beamline provided a pion flux of 1.9- 
2.5 xl0 6 s _1 with e/fj, contamination of 28-58% for mo- 
menta 104-143 MeV/c. The spot size was about 3x3 cm 2 
(FWHM) and the momentum bite was about 6 MeV/c 
(FWHM). Note that the primary proton beam comprised 
pulses of duration 2-4 ns and separation 43 ns, and there- 
fore the arrival of n's, /x's and e's were time-separated. 

Beam particles were registered in two plastic scintilla- 
tors before traversing the target. The upstream counter 
(Bl) was 6.0 cm diameter and 0.3 cm thick and the down- 
stream counter (B2) was 6.0 cm diameter and 0.15 cm 



thick. The targets were disks of polyethylene (CH2) and 
carbon (C) of diameter 7.0 cm and thickness 1.2 cm. 

Neutral pions were identified by detecting the two pho- 
tons from ir° — ► 77 decay using the RMC spectrometer. 
The CH2/C targets were located about 40 cm upstream 
of the geometrical center of the RMC spectrometer, thus 
permitting the detection of forward-going photon pairs 
from forward angle CEX. The photons were detected by 
7 — ► e + e~ conversion in a cylindrical lead converter of 
thickness 1.7 mm, radius 24.7 cm and length 60.0 cm. 
The resulting e + e~ pairs were tracked in a low mass 
cylindrical multi-wire chamber and a large volume cylin- 
drical drift chamber. A four- fold segmented, doublcd- 
layered ring of veto scintillators (the A/ A' ring), located 
within the Pb radius, was used to reject backgrounds 
from charged particles. A sixtcen-fold segmented, single- 
layered ring of trigger scintillators (the D ring), located 
outside the drift chamber radius, was used to identify 
e s from 7-ray conversions. A large solenoid provided a 
2.206 kG axial magnetic field thus enabling the momen- 
tum analysis of charged particles. 

The 7r° trigger required coincident hits in the two beam 
counters, zero hits in the A/A'-ring, and >3 hits in the 
D-ring. The tt° trigger also imposed a requirement of > 3 
cell hits in drift chamber layer 2 and > 6 cell clusters in 
drift chamber layers 3+4. In addition, a beam trigger was 
derived by prescaling the B1-B2 coincidences and used to 
continuously monitor the beam properties. On fulfilling 
either trigger the acquisition read out the chamber hits 
and times and the scintillator ADCs and TDCs. 

A stopping tt~ beam and liquid H2 target were used 
for the spectrometer calibration via at-rest Tt~p — > 7r°n 
using the well-known Panofsky ratio [15]. The stopping 
beam had a momentum of 81 MeV/c and the hydrogen 
target was a cylindrical cell of 16 cm diameter and 15 cm 
length with 0.25 mm thick Au walls. 

Our data was collected at six mementa with central 
values 103.8, 112.9, 118.3, 123.8, 134.4 and 142.8 McV/c. 
At each momenta setting we made measurements with 
the CH2, C transmission targets and an empty target 
setup, and collected all scattering angles simultaneously 

The analysis was organized as follows. First, we ap- 
plied beam, tracking, photon and tt° cuts to identify the 
photon pairs from ir° — > 77 decays. Next, for each n~ 
momentum setting and ir° angle range, the number of 
1 H(ir~,ir°) events was extracted via a fitting procedure 
using the CH2, C, and empty target momentum spectra. 
Finally, from a simulation of the in-flight CEX detection 
efficiency, and the measurement of the at-rest CEX de- 
tection efficiency, we normalized the 1 H(tt~ ,tt°) events 
and determined the cross sections. 

The event reconstruction first assembled chamber hits 
into candidate tracks, then paired e + -e~ tracks into can- 
didate photons, and finally paired photons into candidate 
7r° — > 77 events. A beam cut imposed requirements on 
the minimum amplitudes of the pulse heights in the two 



beam counters. A tracking cut imposed requirements on 
the minimum points and the maximum variances in the 
fits to the helical e + , eT trajectories in the tracking cham- 
bers. A photon cut required that the e + and e~ tracks 
converge at the Pb converter, and a tt° cut required that 
the photon-pair converge at the target. 

The data at each momentum were divided into four an- 
gle bins with scattering angles 0-10, 10-20, 20-30 and 30- 
40 degrees, thus yielding twenty-four data points. Rep- 
resentative 7r° momenta spectra for the CH2 and C tar- 
gets are shown in Fig. 1. The CH2 spectra show a peak 
from the l Yi{it~ ,it°) reaction and a continuum from the 
12 C(7r - ,7r°) reaction. At forward angles the 1 H peak 
and 12 C background were kinematically separated while 
at larger angles the contributions were overlapping. 

To obtain the number of H(7r~ , ir°) events we per- 
formed fits of the measured CH2 spectra to the sum of 
a simulation of the 1 H(7r _ ,7r°) lineshape and the mea- 
surement of the (7r~,7r°) background. Our benchmark 
fits (e.g. Fig. 1) involved one free parameter, the nor- 
malization constant for the 1 H(ir~ , it°) lineshape. The 
background spectra for each angle-momentum bin was 
fixed by combining carbon and empty target spectra in 
accord with the target ir~ exposures and areal densities. 

Additional fits were made to quantify the systematic 
uncertainties associated with the l Y[(it~ ,it°) lineshape, 
(7r~,7r°) background, etc. In one study we allowed the 
position of the *H peak to be varied in the fit, and found 
a mean deviation of the *H yields from the benchmark 
fits of 0.51 a (standard deviations). In a second study 
we allowed the amplitude of the background spectra to 
be varied in the fit, and found a mean deviation of the 
*!! yields from the benchmark fits of 0.37 a. In a third 
study we compared the results of fits either including 
or excluding a small correction for the 1.5-2.5 MeV/c 
differences in the CH2-C target energy loss, and found a 
mean deviation of the *H yields of 0.22 a. For the final 
uncertainties in the H yields we took the quadrature sum 
of the statistical errors from the benchmark fits with the 
maximum deviations from the supplementary fits. 

The number of incident pions was determined by 
counting the B1-B2 coincidences and measuring the beam 
composition. The beam composition was determined us- 
ing the beam trigger data and was based on the differ- 
ences in the flight times and the pulse heights of the it's, 
/z's and e's. The largest uncertainty in the 7r~ flux de- 
termination was the beam composition measurement. 

To derive the l H(n~ , ir°) differential cross sections 
we required the spectrometer response function. The 
response function determines the probability of recon- 
structing a neutral pion of true momentum px and angle 
9t with measured momentum pm and angle 6m ■ It was 
computed via a simulation using the GEANT3 code [16]. 
The simulation incorporated the full geometry of the ex- 
perimental setup and included the detailed interactions 
of the outgoing photons, conversion electrons, etc. The 



same analysis package was used to process both simu- 
lated data and measured data. 

The absolute normalization of the response function 
was obtained by comparing the measurement and simu- 
lation of at-rest CEX for the stopping ir~ beam and the 
liquid H 2 target. The energy-angle distributions of 7- 
pairs from at-rest CEX obtained from the measurement 
and the simulation were in excellent agreement. These 
comparisons were made for different trigger definitions, 
cut parameters, etc., in order to fully test the detector 
simulation. However, the absolute efficiency for tt° detec- 
tion was found to be smaller in the measurement than the 
simulation by a factor, F = 0.85 ± 0.03. This difference 
is suspected to arise from contributions such as chamber 
noise and multiple pulsing that were absent in the sim- 
ulation, and was observed in previous experiments with 
the RMC spectrometer. Importantly, at the level of the 
quoted uncertainty ±0.03, the factor F was independent 
of the trigger condition, cut parameters, photon energies 
and target position. 

During production running the relative acceptance 
was monitored by measurements of at-rest CEX using 
a 2.5 cm thick CH2 target and an 81 MeV/c stopping 
■k" beam. The acceptance was found to be constant to 
about ±5% over the duration of the experiment. Given 
the uncertainties in the determination of the absolute ac- 
ceptance and the monitoring of the relative acceptance, 
an uncertainty of ±10% was conservatively assigned to 
the overall normalization of the cross sections. 

The measured cross sections are averages over the 
momentum dispersion of the beam and the momentum 
spread in the target. The beam momentum dispersion 
was obtained by analysis of the time-of-flight spectra 
of the incident beam particles. The target momentum 
spread was obtained by a GEANT3 simulation of the 
beam interactions in the target material. The typical 
beam dispersion was 2.8 MeV/c (a) and the typical tar- 
get spread was 2.6 MeV/c (a). A correction (of 0.1-0.4 a 
in magnitude) was applied to convert the bin-averaged 
cross sections to bin-centered cross sections. 

Our final results for the center-of-mass (CM) differ- 
ential cross sections versus CM scattering angle and w~ 
laboratory momentum are plotted in Fig. 2. The error 
bars include the combined uncertainties in the tt~ flux, 
1 H(7r~,7r°) yields, angular resolution and momentum in- 
terval, but omit the ±10% normalization uncertainty. 

Our data and earlier data near the minimum are in dis- 
agreement. Compared to Fitzgerald et al. [2] our cross 
sections are roughly factors of two smaller at momenta 
<120 MeV/c. This disagreement is beyond the quoted 
experimental uncertainties and - as discussed in Rcfs. 
[3,13] - may be related to the normalization procedure 
and the background subtraction in the Fitzgerald et al. 
experiment. At 113 MeV/c our results are also ~30% 
smaller than the preliminary results of the Isenhower et 
al. experiment [3], although this discrepancy is statisti- 



cally less significant. 

Also plotted in Fig. 2 are the predicted cross sections 
from Gibbs et al. [4,5], Matsinos et al. [7] and the GWU 
group [12]. The Gibbs et al. analysis used a potential 
model where the parameters of the 7rN interaction were 
obtained by fitting the 7r ± p elastic dataset only. The 
predictions by Matsinos et al. and the GWU group were 
derived from partial wave analyses of the 1 p elastic 
dataset. We stress the GWU results in Fig. 2 are a 'spe- 
cial solution' based on T<100 MeV elastic data only - 
and differ slightly from their well-known SP06 solution 
[12]. All analyses include the isospin violation arising 
from Coulomb effects and hadronic mass splittings. 

We first discuss the comparison among the Gibbs et 
al. and Matsinos et al. predictions and the CEX data - a 
comparison [4-7] that originally motivated claims of large 
isospin breaking effects in 7rN interactions. Fig. 2 shows 
our CEX data and these predictions are in quite good 
agreement, with Matsinos et al. being in agreement at all 
six momenta and Gibbs et al. being in agreement at the 
four lower momentum values and differing from our data 
by about 10-20% at the two higher momentum values. 
In contrast, the original comparisons of these predictions 
to the Fitzgerald et al. data showed large discrepancies. 
We therefore conclude our data contradict the original 
evidence for unexpected isospin violating effects. 

However, the GWU results in Fig. 2 suggest some 
caution is necessary. Their results differ markedly from 
Gibbs et al. and Matsinos et al. in the dip region. Note 
these analyses do differ in details; for example in the 
omission and the normalizations of certain n p datasets. 
As mentioned by Matsinos et al. [7], such differences 
could result in small discrepancies in the s-and p-wave 
partial amplitudes that develop into large discrepancies 
in the CEX cross section. Thus final conclusions on 
isospin breaking will likely require a better understanding 
of such discrepancies. 

Fig. 3 shows the extrapolated zero-degree cross sec- 
tions derived from our data and Fitgerald's data. Al- 
though the extrapolated cross sections have larger un- 
certainties, they do allow for the direct comparision of 
the two experiments. Like Fitzgerald et al, our extrap- 
olations have assumed that the cross sections could be 
parameterized as second-order polynomials in cos 0cm- 
Fig. 3 shows our data and Fitzgerald's data are in clear 
disagreement for momenta below 120 MeV/c but closer 
agreement for momenta above 120 MeV/c. Moreover, 
while our data are reasonably consistent with Gibbs et 
al. and Matsinos et al. at all momenta, the Fitzgerald et 
al. data are obviously inconsistent with their predictions 
below 120 MeV/c. The GWU results are in greatest dis- 
agreement with our data at momenta <120 MeV/c and 
Fitzgerald's data at momenta >120 MeV/c. 

Last, we performed a fit of our zero-degree extrapola- 
tions to a second-order polynomial in the 7r~ momentum 
and found the minimum to be T m i n — 41.9 ± 0.9 MeV. 



Our result is in reasonable agreement with Gibbs et 
al. and Matsinos et al. which yield T m ; in ~ 43 MeV 
but in significant disagreement with the GWU result 
T mm ~ 47 MeV. The Fitzgerald et al. value T min = 
45.1 ± 0.5 MeV and our value T min = 41.9 ± 0.9 MeV are 
in marginal disagreement. 

In summary, we have measured the CEX differential 
cross section at six momenta (104-143 MeV/c) and four 
angles (0-40 deg.) using the TRIUMF RMC spectrome- 
ter. At momenta below 120 MeV/c our cross section mea- 
surements are considerably smaller than the published 
data of Fitzgerald et al. and somewhat smaller than the 
preliminary data of Isenhower et al. Interestingly, our re- 
sults are consistent with the cross sections derived from 
the low energy 7r ± p elastic dataset by Gibbs et al. and 
Matsinos et al., and - unlike the original comparisons be- 
tween these analyses and earlier data - show no evidence 
for unexpected isospin breaking effects. Further investi- 
gations to understand the differences between the CEX 
cross section predictions would be valuable. 
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FIG. 1. The 7r° momentum spectrum for an incident n~ 
momentum 104 MeV/c and scattering angles 10-20 deg. The 
data points (filled circles) are the measured CH2 spectrum 
and the solid line is the fit to the sum of the measured (ir~ ,n°) 
background (light gray fill color) and the 1 H(tt~ ,tt°) lineshape 
(dark gray fill color). 
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FIG. 2. The ix~-p — » -K°n CM differential cross section ver- 
sus central momentum and CM angle. The data points are 
our results. The solid line is the prediction of Gibbs et al. 
[4,5], the dashed line is the prediction of Matsinos et al [7], 
and the dotted line is the prediction of the GWU group [12]. 
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FIG. 3. The zero deg. n~p — ► 7r°n CM difTerential cross 
section versus laboratory momentum. The filled circles are 
our data, the open circles are Fitzgerald's data. The solid 
line is the prediction of Gibbs et al. [4,5], the dashed line is 
the prediction of Matsinos et al. [7], and the dotted line is the 
prediction of the GWU group [12]. 



